Spin-precession-assisted supercurrent in a superconducting quantum point contact 

coupled to a single-molecule magnet 
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The supercurrent of a quantum point contact coupled to a nanomagnet strongly depends on the 
dynamics of the nanomagnet's spin. We employ a fully microscopic model to calculate the transport 
properties of a junction coupled to a spin whose dynamics is modeled as Larmor precession brought 
about by an external magnetic field and find that the dynamics affects the charge and spin currents 
by inducing transitions between the continuum states below the superconducting gap edge and the 
Andreev levels. This redistribution of the quasiparticles leads to a non-equilibrium population of 
the Andreev levels and an enhancement of the supercurrent which is visible as a modified current- 
phase relation as well as a non-monotonous critical current as function of temperature. The non- 
monotonous behavior is accompanied by a corresponding change in spin-transfer torques acting 
on the precessing spin and leads to the possibility of using temperature as a means to tune the 
back-action on the spin. 

PACS numbers: 75.76. +j, 74.50. +r, 75.50.Xx, 75.78.-n 



I. INTRODUCTION 

Spintronics devices, which utilize the spin degree of 
freedom, have already revolutionized read-out technol- 
ogy used in hard drivesi^ In conventional spintronics, 
the transport properties of a device typically depend on 
the relative orientation of the spins of electrons with 
respect to a reference, which may be a magnetic field 
or a magnetization direction of a ferromagnetic layer* 3 - 
The challenge of downsizing electronic devices has lead 
to the study of transport properties of non-magnetic 
single-molecule devices such as diodes^ and transistors^— 
as well as devices containing single-molecule magnets 
(SMMs)4 The interest in SMMs stems from their long 
relaxation times at low temperatures 9 and their display 
of a wide range of quantum physics phenomena^ 9 ' 10 Stud- 
ies on SMM devices include for instance three-terminal 
devices?^ - — supramolecular spin valves^ and inelastic 
tunneling spectroscopy^ 

Interesting spin phenomena may occur when ferro- 
magnets are combined with superconductors (see and 
references therein). Cooper pairs in a conventional su- 
perconductor have spin-singlet pairing which, if the su- 
perconductor is interfaced with a ferromagnet, extend 
into the ferromagnet. However, the exchange field in- 
side the ferromagnet tries to align the two spins of the 
Cooper pairs and hence breaks the Cooper pairs apart 
resulting in a rapid decay of the superconducting cor- 
relations inside the ferromagnet. For the same reasons, 
the critical current of a Josephson junction with a fer- 
romagnetic layer sandwiched between the two supercon- 
ductors decays rapidly with increasing thickness of the 
ferromagnetic layeri2i~— On the other hand, if weakly 
ferromagnetic interfaces with magnetization directions 



differing from the magnetization direction of the fer- 
romagnetic layer are inserted, the spin-singlet correla- 
tions may be transformed into spin-triplet correlations 
which can survive over a long range within the ferro- 
magnet layer * 25-29 As a result of this non-collinear mag- 
netization of the ferromagnetic layer, the critical cur- 
rent decays similarly to a supercurrent in a non-magnetic 
metal with increasing junction length] 30 ' 31 Also interac- 
tion between spin and charge supercurrents has attracted 
attention^ 2 - - — Recently, interest in coupling between the 
dynamics of magnetic moments and Josephson currents 
has increased^ - — 

In Refs. and |H, it was found that spin-singlet 
to spin-triplet conversion can be generated by a nano- 
magnet such as a SMM or a ferromagnetic nanoparti- 
cle coupled to a Josephson junction^ provided that the 
magnetization direction of the nanomagnet processes^ 5 - 
The spin-triplet correlations enable spin currents to exist 
close to the junction interface despite the s-wave nature 
of the superconducting leads. The spin currents generate 
a spin-transfer torque acting on the nanomagnet and its 
effect may be measurable in a ferromagnetic resonance 
(FMR) exp eriment^ as was suggested in Ref. Refs. 
142 and [43| focused on charge and spin currents as well 
as spin-triplet correlations at zero temperature and in 
equilibrium junctions. Here, we instead investigate the 
effects of the spin precession on the critical current and 
find that the critical current is enhanced at high tem- 
peratures. This enhancement is due to a redistribution 
of the population of states caused by the precessing spin 
and is also visible in the current-phase relation (CPR). 
The redistribution of quasiparticles also affects the spin 
currents, which change abruptly as a function of temper- 
ature when the critical current is enhanced. 
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FIG. 1. (Color online) (a) Two superconducting leads are coupled over the spin of a nanomagnet. The hopping amplitude 
for spin-independent tunneling is given by Vo, while Vs is the coupling between the tunneling electrons and the spin of the 
nanomagnet. (b) The spin of the nanomagnet, S, precesses with the Larmor frequency, ujl, due to an external effective magnetic 
field, H , that is applied at an angle •& with respect to the orientation of the spin. The density of states in a rotating frame 
(see definition in the text) is plotted for a transparent junction (vt = 1) with hopping strengths (c) vo = vt cos(0.l7r/2) ss 
0.9877, v s = wt sin (0.l7r/2) sa 0.1564 and (d) v = «tcos(tt/2) = 0,v s = v T sxa(ir/2) = 1 at temperature T/T c -> 0. The 
Andreev levels, which can be seen as sharp states inside the gap, have here for visibility reasons been given an artificial 
broadening, + , and the continuum density of states is normalized to 1 away from the gap edge. The upper (lower) Andreev 
levels, e + (e~), are given an effective Zeemann splitting (e + ' - ' — > e^^~ ) by the spin precession. The spin precession couples 
the Andreev levels and the continuum states as well as the states e^'~' and ( e ) The current-phase relation (CPR) is 

plotted for Vq = cos(0.l7r/2), v$ = sin(0.l7r/2) at temperatures (solid lines from top to bottom) T/T c = 7.5 • 10~ 4 , 0.25, 0.50, 
0.75 and 0.85. The dashed line shows the CPR for (v = l,v s = 0) at temperature T/T c -¥ 0. (f) The CPR for a junction with 
(v =0,v s = 1) at temperatures (from bottom to top) T/T c = 7.5 • 10~ 4 , 0.25, 0.50, 0.75 and 0.85. In both (e) and (f), the 
CPRs develop a peak as the temperature increases. Figures (c) - (f) are calculated for precession angle $ = 7r/4 and frequency 
ljl = 0.2A. 



The paper is organized as follows: In section [TT1 the 
model of the superconducting point contact containing 
the molecular magnet is described. The approach to solv- 
ing the spin-active boundary condition created by the 
molecular magnet is then described in section IIII1 Sub- 
sequently, the results are discussed in section IIVI start- 
ing with the Andreev levels and their population. Then, 
the charge currents and spin currents are discussed along 
with the induced spin-triplet correlations associated with 
the spin currents. Finally, the results are summarized in 
section [V] Appendix [X] details some of the calculations 
in section ITVl 



II. MODEL 

Here, we consider the same junction model as in Refs. 
l42l and |H, i.e. two superconducting leads coupled over 
a nanomagnet, schematically depicted in Fig. HJa). The 
nanomagnet may be a SMM or a ferromagnetic nanopar- 
ticle and is subjected to an external effective magnetic 
field, H , which couples to the spin of the nanomagnet, 
S, through the Hamiltonian T-Lb — 7 H ■ S, where 7 



is the gyromagnetic ratio (see Fig. [IJb)). We treat the 
spin classically and assume that the effective magnetic 
field is applied at an angle i? with respect to the spin. In 
this classical limit, the effective magnetic field results in 
a torque acting on the spin, S = —jH x S, causing the 
spin to precess with the Larmor frequency, u>l = l\H\. 
The effective magnetic field is taken to be applied along 
the z axis and hence the spin, S, can be parameterized as 
S(t) = Ses(t) where S = \S\ is the magnitude of the spin 
and es(t) = cos(cji^) sint?e x +sm.(u)Lt) sin + cos -de z . 
The two superconducting leads are then coupled over the 
nanomagnet leading to a complete Hamiltonian given by 
H = Hieads + H B + Hrr^ where Hi ea ds describing 
the leads is the BCS Hamiltonian with the superconduct- 
ing order parameter A(T)e ±lv ^ 2 , where "+" ("— ") refers 
the right (left) lead. "Ht is the tunneling Hamiltonian, 
U T = J2ka;k'a> c L,ka V ka-k' a' CR,k'a' + h.c , where the hop- 
ping amplitude is given by Vka-.k'a' = (Vb<W + V s (S(t) ■ 
c)aa')S{k — k') and er — [a x ,a v , <r z ) consists of the Pauli 
spin matrices. The first term of the hopping amplitude 
describes spin-independent tunneling, while the second 
term corresponds to spin-dependent tunneling and de- 
pends on the instantaneous direction of the spin, S. The 
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projection of S along the z axis leads to different tunnel 
probabilities for spin-up and -down quasiparticles while 
the transverse components of the spin lead to spin flips. 
Moreover, we assume that the tunneling quasiparticles 
are affected only by the exchange interaction with the 
nanomagnet and that the external magnetic field does 
not couple to the leads or directly to the tunneling quasi- 
particles. 



III. APPROACH 

We use non-equilibrium Green's functions in the qua- 
siclassical approximation 4 ^— to calculate the current 
through the junction ! 52 ' 53 The quasiclassical Green's 
functions are propagators for quasiparticle on classical 
trajectories. An incoming propagator describes quasi- 
particles with trajectories leading into the junction in- 
terface, i.e. vf ■ n < where n is the surface normal and 
vp is the group velocity of the electron-like quasiparti- 
cles. Conversely, outgoing propagators describe quasi- 
particles having trajectories leading out from the inter- 
face, vp-n > 0. The spin-singlet superconductivity of the 
uncoupled leads can be described by spin-scalar Green's 
functions. However, to accommodate for the effects of 
the spin-dependent scattering caused by the nanomag- 
net, the Green's functions are parametrized as consisting 
of both spin-scalar (s) and spin-vector (t) components 



c ,54 



9 



X 



X 



if ' a (f* + ft 
ft • cr ) 9e -Py(9t • <T ) cr y 



(1) 



where X stands for retarded (R), advanced (A) and 
Keldysh (K) Green's functions and the hat (" "") denotes 
matrices in Nambu-spin space. The components and 
f t are the spin-singlet and spin-triplet components of 
the anomalous Green's functions, respectively (similarly 
for gf and gf). 

The precessing magnetization of the nanomagnet con- 
stitutes a time-dependent, spin-active boundary condi- 
tion that can be solved by applying the unitary transfor- 
mation 



W(t) = 











-t<T z 



(2) 



which results in a transformation to a rotating frame of 
reference where the Fermi surfaces of the spin-up and 
-down bands arc shifted with This transforma- 

tion also shifts the gap edges as follows: The upper gap 
edge, A+(T) = A(T), is shifted as A+(T) -> A+ 4 (T), 
where A+(T) = A(T) - auj L /2 and a = 1(-1) for 
spin | U)- The lower gap edge, A~(T) = -A(T), is 
correspondingly modified as A~(T) — > A^(T), where 

A~(T) = -A(T) - auj L /2. In this rotating frame, the 
precessing spin, S(t), now appears static and points along 
the direction eg = cosi?e z +sm$e x . Replacing the hop- 
ping amplitudes Vq and Vs by their Fermi surface limits, 



v = ttNfVq and vs — nNpSVs where Np is the nor- 
mal density of states at the Fermi energy, the hopping 
element in Nambu-spin space reads 



vo + v s e s ■ a 

v 







(3) 



and v — vl in Keldysh space. 

As a first step towards solving the boundary condition, 
the interface is treated as an impenetrable surface and the 
quasiclassical Green's functions are found by solving the 
Eilenberger equationi2r— in each lead separately. The 
resulting solutions, denoted by g° with a = L,R indicat- 
ing the left or right lead, are then connected across the 
interface using a quasiclassical t-matrix equation^ 5 - - — 



i a (e) = f Q (e) + f Q (e)5°(e)t a (e), 



(4) 



which takes the hopping Hamiltonian Ht into account 
through 



\/b{z) = vg R/L (e)v 



(5) 



In this way, the full propagators for the junction can be 
obtained as 

= 9° a (e) + (9° a (e) ± ™i)i a {e){f a {z) T ml), (6) 

where g 1 ^ denotes the incoming (outgoing) propagator. 
The difference between the incoming and the outgoing 
propagators then give expressions for the charge Joseph- 
son current and the spin Josephson current per conduc- 
tion channel as 

f a h (t) = ^J ^Tr[f 3 (^<M) -s^M))]; (7) 

JaW = \J ^Tr{T 3 <T(ti<{e,t)-g%<(e,t))} (8) 

where T3 = diag(l,— 1) and a = diag(cr, — a y <ja y ). 
The lesser ("<") Green's functions is obtained as g K = 
\{g K — g R + g A )- The details are described in Refs. l42l 
andlH 



IV. RESULTS 

Scattering processes between two superconductors may 
lead to constructive interference and the appearance of 
Andreev levels^ In the presence of a precessing spin, a 
tunneling quasiparticle may gain (lose) energy uil while 
simultaneously flipping its spin from down (up) to up 
(down). These additional tunneling processes lead to a 
modified Andreev level spectrum whose details depend 
on a number of parameters: The ratio vq/vs determines 
whether the junction is in a or a f state, depending 
on if vo/vs > 1 or v /vs < 1, respectivel y 42 ' 43 The Lar- 
mor frequency, wj,, determines the amount of energy a 
quasiparticle may gain or lose during tunneling across 
the junction. The precession angle, determines the 
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amount of scattering between the spin-up and -down 
bands. The population of the Andreev states is modi- 
fied by the temperature, T, but also by the scattering 
processes generated by the precessing spin. This modifi- 
cation of the Andreev level population at finite temper- 
ature is the focus of this paper. 



A. Andreev levels 

In the case of a static spin, for which one can take 
u>l = and $ = 0, the Andreev levels are given by 



±AW1 - D sin 



2 



T>s cos 2 



where 



1 + 2{v 2 + v s ) + {vl 



4) 2 



(9) 



(10) 



are transmission probabilities determined by the hopping 
amplitudes vo and vs- The two signs of e ± describe the 
two Andreev level branches; one that exists below the 
Fermi surface, e~ , and one that exists above the Fermi 
surface, e + . From Eq. ©, it is clear that the junction is 
in the state if vq > vs, and in the ir state if vq < vs^ 
In the limit of vs -+ 0, but with an arbitrary preces- 
sion frequency, the Andreev levels in the rotating frame 
display an effective Zeemann splitting, — > ' 
The two branches are then given by 



't.4- 



A 



1 - V sin" 



(I) 



2A 



-A 



1 - V sin z ( % 



Note that these states are eigenstates only when there is 
no spin-flip scattering and that spin-flip scattering leads 
to scattering between and \ The Andreev 

levels are visible in the density of states as sharp subgap 
states. The time-averaged density of states is evaluated 
as 



Pa(e,<p) 



1 

'8tt 



3{Tr[f 3 (^(e,^+^(e,^))]}, 

(12) 

where g& (e, p) is the diagonal component of 
ga°' R (e,(p,t) and is time independent (see Ref. In 
Fig- QI C ), tne density of states is plotted in the limit 
vs *C Vo- The sharp states inside the superconducting 
gap are the Andreev levels, which are split into and 

ej. These Andreev levels are well described by Eq. pT|) 
and their splitting is given by w^. 

Increasing the spin-dependent hopping amplitude, vs, 
results in scattering of quasiparticles between the split 
Andreev levels belonging to the same branch. In general, 
the Andreev levels can be described by 

= ±Ai/l + $(v ,v s ,uj L ,fl,<p) + Z a (v ,v s ,u} L ,fl,<p), 

(13) 



where an analytical expression for $ can be found as 
*(vo,v s ,ul,0,<p) = -V ($) S m 2 (|) (14) 

and the transmission probabilities, Vs and Vo, depend 
on the precession angle $. Defining 6 = 1 + 2(vq + v s ) + 
(Vq — v s ) 2 , the transmission probabilities are given by 



o i'I , cos(2z?) 2v§sin 2 (2tf) 



Vs(&) 2v s\ § + 6-4v s sm 2 ti ' (6 - Av% sin 2 0)2 



and 



Vo{&) = 2v% Q + 



cos (2i?) 
6 - Av 2 s sin 2 <d 



(15) 
(16) 



The function S CT (?;o, vs, ojl, p) provides the Zeemann 
splitting which in the general case also depends on i9. An 
analytical expression for can in principle be obtained, 
but is too involved to allow for simple analytical anal- 
ysis. The effects of this term are, however, numerically 
analyzed below. The density of states for a junction with 
(vq = 0,vs = 1) is shown in Fig. Hfd). The Zeemann 
splitting in this case, vs 3> vo, can be approximated with 
~ lol cosfl. 



B. Enhancement of the critical current 

Each Andreev level carries a certain amount of cur- 
rent that is weighted by the Andreev level occupation. 
In equilibrium for a static spin, the amount of current 
each Andreev level, e ± , carries is {2e/Ti)de ± / 'dp while 
the population of the quasiparticle states is given by 
the Fermi distribution function, <fto ' 59 ' 60 In this equilib- 
rium situation, the lesser Green's functions in Eq. ([7]) 

can be written as ga^'a (s, p) = —2irp]l iya {e,p)4>o{^) 
using the partial density of states Pu° aa > f) — 

(~s)^{te"°' S )™'( £ '^}' Hence, in equilibrium, the 
charge current is given by 



f*M = s 



(V - Vs)Asm(p) tanh[e+(v?)/2T] 



1 - V sin 2 (<^/2) - V s cos 2 {p/2) 



(17) 



In a non-equilibrium situation where a superconduct- 
ing point contact contains a spin precessing with a finite 
frequency, the current- pha se relation is modified. As was 
shown in Refs. HH and l43l the Josephson charge current 
is time independent and in the case of dominating spin- 
dependent tunneling and zero temperature, the CPR ex- 
hibits abrupt jumps as a function the superconducting 
phase difference, p. As the temperature is increased, the 
abrupt jumps are smoothed out and a new step at a phase 
difference p p develops as can be seen in Fig. [IJe) and (f). 
This step, which consequently gives a peak in the CPR 
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FIG. 2. (Color online) (a) Population of the lower Andreev 
level e^jT (black lines) and the upper Andreev level (red 
lines) as a function of phase difference, ip, for a junction with 
hopping amplitudes vo = cos(0.l7r/2) and vs = sin(0.l7r/2), 
precession frequency ujl = 0.2A(T = 0) and tilt angle •& = 
tt/4. At temperature T/T c = 7.5 ■ 1CT 4 (dotted lines), the 
population of the Andreev levels is unaffected by the phase 
difference. At higher temperatures, T/T c = 0.25, 0.50, 0.75 
and 0.85 (solid lines), the population changes abruptly at a 
phase difference ip p leading to a jump in the CPRs shown 
in Fig. [TJ (b) Density of states for the junction in (a) at 
temperature T/T c = 0.85. The dashed lines mark A ± (T) = 
±A(T). The dotted lines denote A±(T) = A ± (T) - auj L /2. 
The population step in panel (a) occurs at the phase difference 
ifip given by the condition ujl ~ A*(T)| — \sf(T,(p)\. 



at ip p , is the result of an enhanced Josephson current at 
high temperatures for phase differences in the interval 
< ip < ip p for junctions and for ip p < ip < n for 7r 
junctions. To understand this current enhancement, we 
now turn to the population of the Andreev levels. 

The population of the spin band a in a non-equilibrium 
situation can analogously to the equilibrium case be de- 
fined as 

cfl°><(e,ip) = gZ a < {e, V )/[-2^^{e^)]. (18) 

In Fig. H[a), the occupation of the Andreev levels el 
(black lines) and (red lines) in a junction with (vq = 



cos(0.l7r/2), vs = sin(0.l7r/2)) is plotted as a function of 
phase difference, ip, for temperatures T/T c = 7.5 • 10~ 4 
(dotted lines), 0.25, 0.50, 0.75 and 0.85 (solid lines). At 
temperature T — >• 0, the lower Andreev level, e^jT, is fully 
occupied while the upper Andreev level, e^, is unoccu- 
pied. At higher temperatures, the lower Andreev level is 
not fully occupied while the upper correspondingly has 
a finite occupation. In addition, there is also an abrupt 
change in the population corresponding to the jump in 
the CPR at phase difference (p p . The abrupt change in 
the population is an effect of the spin-flip scattering pro- 
cesses in which a quasiparticle interacting with the pre- 
cessing spin may gain or lose energy u) l ■ These scattering 
processes couple the Andreev levels with the continuum 

+(-)/ 



states provided that w L > '(T)\ 



\et { '\, see Fig. 
Hfb). The coupling causes quasiparticles to be promoted 
to the lower Andreev level from the continuum below the 
gap edge. Quasiparticles in the upper Andreev level are 
similarly scattered into the upper continuum. If the lower 
(upper) Andreev level is not completely filled (unoccu- 
pied), which is the case at finite temperature, the cou- 
pling leads to a repopulation (emptying) of the Andreev 
level similarly to the repopulation of Andreev levels due 
to microwave radiationi 61 i 62 This process enhances the 
supercurrent as is shown in Fig. [TJe) and (f). The peak 
in the CPR at phase difference ip p can hence be found 

from the condition ujl — (T)\ — \e^ '\, which can 

be obtained from Eq. (ITT]) in the case of vs <C Vo, and 
leads to 



2 arcsin ■ 



UJL 



A)A(T) 



2 A(T))/» 



(19) 



where u>l < £>oA(T)/2. As a consequence, tp p increases 
as the temperature increases and the superconducting 
gap, A(T), closes. 

The enhancement of the charge current for certain 
phase differences may lead to an enhanced critical cur- 
rent if the enhancement is large enough. Fig. [3] 
shows the enhancement of the critical charge current, 
Ij^i^^l/lj^i^L = 0)|, as function of temperature, 
T/T C) for precession frequencies ljl = OTA, 0.2A and 
0.3A. Dashed lines denote hopping amplitudes (v — 
cos(0.l7r/2), vs = sin(0.l7r/2)) while solid lines denote 
(vq = 0, Vs = 1). The enhancement of the critical current 
may be understood as an effective lowering of the temper- 
ature due to the repopulation of the lower Andreev levels 
and the emptying of the upper Andreev levels - the An- 
dreev level population generated by the precessing spin at 

a certain temperature, 4> 1 J°' < {u)l > 0,Xi), corresponds to 
an Andreev level population at a lower temperature but 



with zero precession frequency, 



i/o,< 



(u L = 0,T 2 <T 1 ). 



The critical charge current as a function of preces- 
sion frequency is shown in Fig. [4] for hopping ampli- 
tudes (a) (v — cos(0.l7r/2), vs = sin(0.l7r/2)) and (b) 
(vq = 0, vs = 1). The critical current is enhanced at high 
enough precession frequencies due to the redistribution 
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FIG. 3. (Color online) The spin-precession-assisted charge 
current in Fig. [T] may influence the critical charge current, 
j° . This is the case at high temperatures, as shown here 
where the enhancement of the critical current compared to 
the critical current for a static spin, \j^ h '(u>L,)\/\j2 h '(wl = 
0A)|, is plotted as a function of temperature, T/T c , for 
(vq — cos(0.l7r/2), vs = sin(0.l7r/2)) (dashed lines) and 
(vq = 0, Vs = 1) (solid lines). The precession angle is $ = n/8 
and the temperature behavior of the critical current for a 
static spin is given by jf{uj L = OA) oc A(T) tanh [^P] .— 



of the quasiparticle occupation. At even higher preces- 
sion frequencies, quasiparticles can be scattered between 
the lower Andreev levels, e~, and the upper continuum 
states, similarly to what was found in Refs. and [H. 
These processes lead to a decrease of the critical current 
at precession frequencies ~ 2A(T). 



C. Spin currents and spin-triplet correlations 

In Refs. 42 and 43s, it was found that the spin- 
dependent scattering across the junction produces a spin 
structure of the full propagators g]f . In general, the 
Keldysh-Nambu-spin matrices of the boundary condi- 
tion problem have a spin structure that enables them 
be divided into generalized diagonal matrices, A , spin- 
raising matrices, at, and spin-lowering matrices, X^. 
This spin structure of the g l ^° matrices, in particular, 
is what leads to the non-zero spin current of Eq. ([8]l. 
The spin current exists close to the junction interface 
and decays on the scale of the superconducting coherence 
length. In Refs. and El the spin current was found 
to have a precessing polarization and a term due to the 
spin-dependent Andreev scattering that can be expressed 
as j s a H (^H) x S, which is finite only for temperatures 
T < T c . This term generates a feed-back effect on the 
precessing spin in the form of a spin-transfer torqu o 63 > 64 
given by t h = (j s L H - jJ^)( 7 H) x S^. This spin- 




(a) V A T=0 0>) VN-O 

FIG. 4. (Color online) Critical current as function of pre- 
cession frequency, ujl, for (a) (vq = cos(0.l7r/2), vs = 
sin(0.l7r/2)) and (b) (vq = 0, us = 1), where the preces- 
sion angle is taken to be $ = n/8 and the temperature is 
T/T c = 7.5 -10" 4 (black), 0.25 (red), 0.50 (green), 0.75 (blue) 
and 0.85 (violet). The increase of the critical current at low 
frequencies is an effect of the spin-precession enhancement of 
the current which occurs at frequencies high enough to cou- 
ple the lower (upper) Andreev levels with continuum states 
below (above) the gap edge, i.e. ujl > [Aj(T)| — \e„(T,ip)\. 
The decrease of the critical current at frequencies ujl ~ 2 A 
is due to coupling between the lower (upper) Andreev levels 
with the upper (lower) continuum states. 

transfer torque acts as an effective magnetic field and 
shifts the precession frequency as u L — > + 2j R ], 

where f H = j s L H = -j s RM - 

The superconducting phase difference ip c correspond- 
ing to the critical current, i.e. jf 1 = j ch {f c ), is plotted in 
Fig. [5] for (a) (v = cos(0.1tt/2), v s = sm(0.l7r/2)) and 
for (b) (uq = 0,^5 = 1). As can be seen in the figure, 
the phase ip c changes abruptly as tp c coincides with tp p . 
The abrupt change in ip c is accompanied with an equally 
abrupt change in the spin-current component j s H , which 
is shown in panels (c) and (d) of Fig. [S]for the two choices 
of hopping amplitudes. 

The connection between the spin currents and the 
anomalous Green's functions generated by the precess- 
ing spin can be described by relatively simple relations 
in the case of a small tilt angle, d = 5d. If the t matrix 
and the incoming (outgoing) propagator in the case of 

$ = are t° and §a 1 ^ , respectively, then one can show 
that (see Appendix [A| a small tilt angle 5$ > gives a 
contribution, 5i and 6g, to first order in <5i9 such that 

i a = i° a + St a (20) 

and gt=& i,0 + S~9t- 

The first-order contribution to the t matrix has the spin 
structure 5t a = 5i\ + St^ and hence the contribution to 

the propagators takes the similar form 5ga° — 5ga' l ^° + 

$ga' 1 where 

6§W° = (S a ±mi)5fl a ($ a Ti*l) (21) 
and 5g^/° = (g° a ± ^1)^(5° T ivrl). 
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FIG. 5. (Color online) In panels (a) and (b), the superconducting phase difference, tp c , corresponding to the critical charge 
current, i.e. jf 1 = j ch ((fc), is plotted as a function of temperature, T/T c , for (a) (vo = cos(0.l7r/2), vs = sin(0.l7r/2)) and 
(b) (vo — 0,vs = 1) for precession frequencies ujl = 0.1A, 0.2A and 0.3A. The abrupt jumps in the phase difference occurs 
when the critical current is given by the spin-precession-enhanced current. The abrupt change in ip c leads to an abrupt change 
in the spin-current component jfj (<p c ) which is plotted for hopping strengths (c) (vo = cos(0.l7r/2), vs = sin(0.l7r/2)) and (d) 
(vo = 0, vs = 1). The spin-current component jjj(ip c ) modifies the precession frequency o;_ l 42 ' 43 and according to panels (c) 
and (d), the magnitude of this frequency shift is drastically modified as a function of temperature. The precession angle is 
$ — it/8 in all panels. 



Using Eqs. © and (p?Tj) . one can write the normal and 

anomalous Green's functions of g l J° — + Sg l J° as 
(see Eq. O) 



(0,i/o,X x i/o,X 

9+ og, 
r i o.X OA o,X 

5g{ g- 



and 



c- pilo.X pQA/o,X 



The normalization condition, 



-2 2} 

9 =-TTl, 



(22) 



{fi^ x +f\'°' x -<r)^y=y_ji /0 ,x j ; f i ,x I m) 



(24) 



can be used to relate the anomalous and the normal 
Green's functions. To first order in <5$, the normalization 
condition gives for the retarded (R) and advanced (A) 



components (of both the incoming and outgoing propa- 
gators) 

Sg R/A = p R /A s f R/A + F R/A S JR./ A (2g) 

5gf A = F R/A 5 ff /A + F R/A 6 ff ,A (26) 

where the coefficients F R j A and F±^ A are given by Eq. 
(IA13[) in the Appendix. The Keldysh components are 

8g? = 8f R F' K + 5f A F[ K - 8f R F R G K - 5f A F A G^ 
+5f?F' A + SffFf (27) 

S 9 f = SffPf + 5f A P' K 5f«F R Qf 5f A F A Q« 
5f«P' A + SfJP' R , (28) 

where the retarded-, advanced-type coefficients, P'±^ A , 



F. 



iR/A p'R/A p'R/A 



± > 



and Keldysh-type coefficients, 



piK F^, P'£ Pg, G± and Q%, are given by Eqs. 
(|A16|1 and (|A17[) in the Appendix. 
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The induced spin-triplet components were quantified 
using d vectors in Ref. |42] as these are used to character- 
ize the order parameters of spin-triplet superconductors. 
However, the appearance of induced spin-triplet correla- 
tions is related to the transport properties of the junction 
as well. The charge current in Eq. is given by the di- 
agonal elements of g l Jt 0,< . To first order in the charge 
current is then given by 

j " {t) = ^n J ^ml'^i^-g^isM, (29) 

i.e. the charge current is identical to its value in the d = 
case. The spin current in Eq. (|8]) is given by the triplet- 
components, which also include the off-diagonal elements 
of the normal Green's functions, and since the z compo- 
nent of Sga° is zero, the spin-up and -down components 
of the spin current can explicitly be written as 

(30) 

where the symmetry relations for the retarded, advanced 
and Keldysh Green's functions were used^ 4 - and the " <" 
components were obtained as x K = \{x K — x R + x A ). 
From Eqs. (J25) - (® and f3"0"|. it is clear that the 
anomalous Green's functions induced by the process- 
ing spin are what determine the spin current. Conversely, 
the spin current vanishes, a (t) = 0, for non-existent 

anomalous Green's functions, i.e. f^'^ = 0. 

V. CONCLUSIONS 

In summary, we have studied the non-equilibrium 
transport properties of a superconducting point contact 
coupled to the precessing spin of a nanomagnet. First, 
we analyzed the Andreev levels and their population as 
a function of temperature. It was found that the preces- 
sion of the nanomagnet modifies the Andreev scattering 
in such a way that it leads to scattering of quasiparticles 
from below the superconducting gap edge into the lower 
Andreev levels. These transitions are similar those oc- 
curring in microwave-irrated superconducting quantum 
point contact s 61 ! 62 and lead to a non-equilibrium popu- 
lation of the Andreev levels and an effective cooling of 
the point contact. 

We have also shown that the effective cooling leads to 
an enhanced supercurrent. If the precession frequency 
of the nanomagnet 's spin is large enough, this spin- 
precession-assisted supercurrent has an enhanced criti- 
cal current. The enhancement of the critical current due 
to the spin precession increases with increasing temper- 
ature. 

The scattering across the junction generated by the 
precessing spin leads to additional non-equilibrium pro- 
cesses: Besides the enhanced supercurrent, we have also 
shown that spin-triplet correlations are created as a re- 
sponse to the scattering processes generated by the spin 



precessing at a finite angle. Moreover, we have shown 
that the induced spin-triplet correlations create a spin 
current. As was shown previously ; 42 ! 43 one component 
of the spin current shifts the precession frequency of the 
spin. Here, this spin-current component was studied for 
the same superconducting phase difference that is asso- 
ciated with the critical current and the magnitude of this 
component was found to be drastically modified as a 
function of temperature. The critical current enhance- 
ment is an effect of the non-equilibrium processes taking 
place inside the junction and a measurement of this en- 
hancement would suggest the existence of induced spin- 
triplet correlations since both phenomena have the same 
origin. 

A natural question concerns the experimental con- 
trol over the junction parameters. Aluminum micro- 
bridges can be used to fabricate few-channel super- 
conducting atomic point contacts 6 -^ whose transmis- 
sion eigenvalues are possible to determine via trans- 
port measurements , 69 i 70 The superconducting gap in 
aluminum is ~ 200 /ieV but can be made smaller in 
an atomic point contact. Typical values for preces- 
sion frequencies are in the range of tenths of GHz for 
ferromagnetic-resonance experiments performed on thin 
ferromagnetic films in contact with superconductors^ 6 - 
With an applied external magnetic field of ~ 60 mT, 
the corresponding gyromagnetic ratio is close to that of 
free electrons, 7 = 2(j,b/^- Assuming similar values for 
the nanomagnet, the precession frequency is in the range 
wl/A ~ 0.07 — 0.3 for a junction with a superconducting 
gap in the range A ~ 20 /xeV— 100 ^eV. 

A direct measurement of the current-phase relation can 
in principle be done for atomic point contacts 71 . How- 
ever, this type of measurement requires the addition of a 
SQUID loop to control the phase difference of the super- 
conducting junction and the magnetic flux through the 
loop might interfere with the magnetic control of the dy- 
namics of the nanomagnet. Therefore, measurements of 
the critical current as a function of temperature should 
be the more practical route to finding experimental sig- 
natures of the coupling between a superconducting junc- 
tion and the dynamics of a nanomagnet and could lead 
to more insight into the interplay between superconduc- 
tivity and ferromagnetism. 
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Appendix A: Calculation of spin-triplet correlations 

In the limit of a small precessing angle, $ = the 
contributions to the t matrices and propagators can be 
calculated as follows: Starting with a junction with zero 
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tilt angle, 1? = 0, the hopping element takes the diagonal 
form 



v + v s <J z 

v a + vsa z 



(Al) 



and the corresponding t-matrix equation is 

^) = f°( £ ) + [f°o^o») (A2) 



where a denotes the left (L) or right (R) lead, c/° is the 
unperturbed Green's function and t° L ^ R = v d o g R , L o v d . 

This equation has the solution i° a = [I - f^g"] -1 ?". 
Notice that the t matrix i° has a diagonal form due to 
the spin structure of r°. The incoming and outgoing 
propagators, g % and g°, given by 

9°/'° = S°a + (9°a ± ™i)Wa T W), (A3) 

have the same diagonal form as t° a and yield a spin cur- 
rent j s a = (see Eq. ©). 

If a small tilt angle, •& = <5i?, is introduced, the hopping 
element is modified into 



The incoming and outgoing propagators are then given 



by g, 

and 



/o -0,i/o 

- go 



Sg l J , where 



t,i/o = /xo 



r**/° _ x?,t,i/° 



(g u a ±i7ri)8tl(g° a T^i) (A10) 



\.,i/o 



and ^/ o = (g0±i7ri)^(^T^i). 

Using the normalization condition, Eq. and the 

explicit forms of the normal and anomalous Green's func- 
tions given by Eqs. and (f2"3")l . respectively, the re- 
tarded and advanced normal Green's functions can be 
expressed in terms of the anomalous ones as 



R/A = pB/A Sf B/A + F R/A 5 ~ f E/A 
R/A = pR/A Sf R/A + F R/A S JR/A 



where 



? R/A 
■ ? R/A 



f 0, R/A 
J± 
p),R/A 
J± 



/£0,R/A 
/£0,R/A 



(AH) 
(A12) 

(A13) 



V + V S (T Z + VsStfcTx 

v + v S (J z + v s 5i9<j x 

(A4) 

where it was used that cos(S'd) f=s 1 and sin(M) s» Sd. 
The hopping element, v, can hence be divided into v = 
ii d + + tr, where v d is given by Eq. (|A1[) and both v< 
and are proportional to 6$: 



til) - ( «s*0o+(-) 








(A5) 



With this hopping element, the change in the t matrix 
can be written as 



t - i° 

l a ~ l a 



St a , 



(A6) 



where t° is the solution of Eq. (|A2[) and 6t a is due to 
the tilt angle, <W. Similarly, 



r - r u 



■ST r 



(A7) 



where f ° = v d g R ^ L v d and the tilt angle gives the 

contribution <5r a . However, the diagonal component of 
St a is 8T d L i R — v^g R / L v^ + v^g R , L v^ w to first or- 



der in 5~&. Then, 6T 



5ft + gr± where Sfl /R 



Here, g± and f± are the diagonal matrix components 
given by g£ = g? ± gf and jf = ff ± ff and (, R ' A = 



0, R/A 
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0,R/A 



The off-diagonal Keldysh components 



8g$ = 5f?FL K + 5ffF[ K - 8f R F R G K - 8ffP A G K 



-5f?F' A + Sff F[ R 



+ 

(A14) 



5gf = &f*P+ + 5~ft p ' K - 5~fiF R Ql - 6ffF A QK 
5 fi p + + <5/f PL R (A15) 

where the retarded-type and advanced-type matrices are 



-,/R/A 



= /r M /Ci ; p±'~ = f± "'"ft* (Ai6) 



y/R/A 



Q,R/A 



F. 



I R/A 



JO, R/A if n'R/A 70, R/A 

j± /Ci;-f± - j± /w 



0,A 



with Ci = .g 4 
exchanging — f-> 



<7_' and C2 is obtained from Ci by 
The Kcldysh-type matrices are 



711K 



G± 
= [F-. 



C±0±/Ci;O± 



[F R 



F_ 



±J ( ^±.- r ± 



P 



/A' 



C±4/C2 
= [P 



(A17) 



i°±;Pf 



P^] 



4 
4 



which means that St l a) cx 6$. Inserting Eqs. and where ^o^ is the occupa tion associated with g° ± and Cj 



i7|) into the t-matrix equation and subtracting Eq. ( IA2 
gives 



o,A 



0,A 

3± • 



5t n 



5T a [l + g'X]+K9> a 



(A8) 



to first order in Sfl. Due to the spin structure of 5t a , the 
t matrix is Si a = Si^ + Si^ where 

^ = [i-iM-W W) P + m (AQ) 



We have thus shown that a small tilt angle, 5$, pro- 
duces spin-flip scattering that combined with Andreev 
scattering processes leads to spinful anomalous Green's 
functions, Sf^^, that in turn generate off-diagonal, spin- 
ful normal Green's functions, Sg^^. Consequently, the 
spin current given by Eq. ([5]) is non-zero as a result of 
the induced spin-triplet anomalous Green's functions. 
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